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ABSTRACT 


This  document  presents  a  model  for  estimating  battery  performance  as  a  function  of  the 
power  level  demanded.  General  values  or  data  specific  to  a  particular  battery  may  be  used  to 
generate  the  required  parameters  for  the  model.  All  required  information  may  be  obtained  from 
a  typical  vendor  data  sheet  that  includes  constant-current  discharge  profiles.  The  total  model  is 
constructed  of  three  separate  components:  (1)  a  model  of  cell  open-circuit  voltage  as  a  function 
of  battery  state-of-charge,  (2)  capacity  de-ratmg  for  constant-current  discharge  per  Peukert’s 
equation,  and  (3)  application  of  Peukert’s  equation  to  variable  current  discharge. 
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1.  INTRODUCTION 


As  an  alternative  to  the  now-common  method  of  battery  system  modeling  by  estimating 
overall  values  of  specific  energy  and  specific  power  on  gravimetric  and  volumetric  bases, 
independent  of  the  rate  of  discharge,  the  model  presented  here  is  recommended  for  estimating 
battery  performance  as  a  function  of  the  power  level  demanded.  General  values  or  data  specific 
to  a  particular  battery  may  be  used  to  generate  the  required  parameters  for  the  model,  shown  in 
table  1 .  All  required  information  may  be  obtained  from  a  typical  vendor  data  sheet  that  includes 
constant-current  discharge  profiles.  Additionally,  development  of  the  model  requires  the 
variables  in  table  2.  The  total  model  is  constructed  of  three  separate  components: 

•  a  model  of  cell  open-circuit  voltage  as  a  function  of  battery  statc-of-charge,1 

•  capacity  de-rating  for  constant-current  discharge  per  Peukert's  equation,  and 

•  application  of  Peukert’s  equation  to  variable  current  discharge.1 


Table  1.  Battery  Mode l  Input  Parameters 


Symbol 

Battery  Parameter 

Efuli. 

Fully-charged  potential 

Eexp 

Potential  at  end  of  exponential  zone 

E\om 

Potential  at  end  of  nominal  zone 

Ecut 

Cutoff  potential 

In 

Specified/measured  discharge  current 

m 

Mass 

pe 

Peukert  effect  exponent 

Qexp 

Capacity  discharged  at  end  of  exponential  zone 

Qnom 

Capacity  discharged  at  end  of  nominal  zone 

Qcut 

Capacity  discharged  at  cutoff 

Rc 

Internal  resistance 

V 

Volume 

Table  2.  Additional  Variables 


Symbol 

Variable 

C 

Capacity  discharged 

Eoc 

Open-circuit  potential 

ED 

Energy  density  (volumetric  basis) 

ES 

Specific  energy  (gravimetric  basis) 

/ 

Actual  discharge  current 

P 

Discharge  power 

t 

Time 

T 

Total  discharge  time 

1 


2.  EFFECT  OF  DISCHARGE  CURRENT  ON  EFFECTIVE  CAPACITY 


2.1  PEUKERT’S  EQUATION 

In  his  study  of  lead-acid  batteries,  Peukert2  observed  a  decrease  in  the  available  capacity  of 
the  battery  with  increasing  discharge  current,  which  could  be  modeled  as 

Ipe  •  t  =  constant, 

where  1  is  the  discharge  current,  /  is  the  time  of  discharge,  and  pe  is  a  coefficient  (“Peukert  effect 
exponent”)  greater  than  one.  Rearranging  Peukert's  law  in  terms  of  capacity  C,  Peukert’s 
equation  may  be  written  to  express  the  capacity  (Ci)  for  a  particular  discharge  current  (7|)  as  a 
function  of  a  known  current  and  capacity  pairing  (such  as  the  nominal  values  frequently  reported 
for  commercial  batteries,  designated  here  as  Iq  and  Co): 

c^-er1-  to 

The  value  of  pe  is  empirically  determined  from  experimental  discharge  data. 


2.2  APPLICATION  OF  PEUKERT’S  EQUATION 

To  model  systems  where  the  requirement  is  a  constant  delivered  power,  not  a  constant 
current  discharge,  equation  (1)  must  be  written  in  such  a  way  as  to  allow  computation  of  C\  from 
a  variable  I\.  If  time  t  is  allowed  to  be  an  independent  variable  against  which  the  various  battery 
parameters  may  be  tracked  (i.e.,  the  battery  capacity  Ci  is  discharged  at  a  rate  of  I\  in  time  t|), 
then 


£i  _  (io\ 
Co  \lj 


pe- 1 


 h-ti  


h-ti 


I 


Effective 


to  ‘  to  I  Effective 't  l 
pe- 1 


C Effective  ~  2(  ^ Ef fcctivc  '  At)  . 


(2) 

(3) 


The  total  capacity  that  may  be  discharged,  per  equation  (3),  may  be  set  equal  to  the  nominal 
capacity  of  the  battery  for  determination  of  the  available  capacity  during  a  variable-current 
discharge.' 


o 
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3.  OPEN-CIRCUIT  POTENTIAL  AS  A  FUNCTION  OF  STATE-OF-CHARGE 


A  simple  model  for  determining  battery  open-circuit  potential  from  the  state  of  charge  was 
developed  in  Tremblay  et  al.1  This  model  is  meant  to  simulate  a  battery  as  it  is  charged  and 
discharged;  i.e.,  it  is  applicable  to  secondary  battery  systems.  The  model  specifically  omits 
voltage  recovery  (an  increase  in  voltage  independent  of  the  state-of-charge)  and  assumes  the 
following: 

•  constant  internal  resistance, 

•  discharge  characteristics  of  the  battery  represent  the  reverse  of  charge  (reversible 
process), 

•  no  effect  of  current  on  capacity  (no  Peukert  effect), 

•  no  temperature  effects, 

•  no  self-discharge, 

•  no  memory  effects. 

The  model,  as  presented  by  Tremblay  et  al.,  is 

E°c  =  E<>-  +  (A  •  exP(-fi  •  C))  ' 

where 

^  =  Efull  ~  Eexp  > 

K  _  (Efull  -  g)voM  +  ><(exp(-g  •  Qnqm)  -  l))  •  ( Qcut  ~  Qnom ) 

Qnom 

E0  =  EFuii  +  K  +  (/?c  •  Iq)  —  A  . 


The  model  parameters  are  found  from  published  manufacturer  data  and  by  inspection  of 
constant-current  discharge  curves.  Determination  of  the  several  values  of  E  and  Q  from  a 
discharge  curve  may  be  made  according  to  figure  1 .  Emu  (and  Qmu>  which  is  by  definition 
zero)  would  be  read  from  the  point  labeled  “Fully  charged  voltage.”  Similarly,  Eexp  and  Qexp 
would  be  read  from  “End  of  exponential  zone,”  Esom  and  Qnom  would  be  read  from  "End  of 
nominal  zone,”  and  Ecur  and  Qcut  would  be  read  from  the  end  of  the  curve.  In  the  case  of  an 
abscissa  marked  in  time  and  not  capacity,  as  in  figure  1,  determination  of  the  values  of  Q  is  the 
product  of  the  value  of  the  discharge  current  and  the  time  indicated  at  each  point.  Additional 
details  may  be  found  in  Tremblay  et  al.1 
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Figure  1.  Identification  of  Discharge  Profile  Characteristic  Values1 


Nominal  Current  Discharge  Characteristic 


The  model  was  validated  against  discharge  curves  for  an  individual  Panasonic  cell  of  each 
of  the  following  chemistries  (all  secondary  cells): 

•  lead-acid, 

•  nickel-cadmium, 

•  lithium,  and 

•  nickel-metal  hydride. 


The  authors  of  the  above  model  presented  a  revised  version  of  equation  (4)  in  Tremblay  and 
Dessaint,4  including  equations  tailored  specifically  to  each  of  the  four  secondary  cell  chemistries. 
The  fine  degree  of  increased  accuracy  provided  by  this  improved  model  is  insufficient  to 
outweigh  the  costs  of  sacrificing  a  general  equation  of  cell  voltage  for  chemistry-specific 
formats;  consequently,  this  revised  model  has  been  disregarded  here. 


4.  A  CONSTANT-POWER,  CONSTANT-TEMPERATURE  DISCHARGE  MODEL 

4.1  ORDINARY  DIFFERENTIAL  EQUATION  FOR  CAPACITY 

While  many  battery  cells  may  demonstrate  nearly  full  capacity  discharge  at  currents  near  the 
published  or  nominal  rates,  it  is  desirable  to  make  theoretical  predictions  of  battery  behavior 
when  systems  are  stressed.  Integration  of  a  capacity  de-rating  model2’3  into  a  model  for  potential 
makes  possible  a  more  conservative  evaluation  of  cell  discharge  at  increased  currents.  Further, 
the  model  may  be  rearranged  to  define  a  constant-power  state: 

WpEfJecttve)  _  [b-Vft2  -4-a-P|Pe  (f. 

At  “  [  2  a  J  '  * 


4 


where 


(2(pe-l)\ 

'o  PC  )• 

b  =  [E°~  {qcut- effective)  +  ^  ‘  exp(_£?  '  Effective))]  ’  ('o  ^  )  - 


^Effective]t=0  —  0* 

The  derivation  of  equation  (5)  is  given  in  appendix  A. 


4.2  APPLICATION  TO  SIZING  ESTIMATIONS 

Integration  from  zero  time  until  C Effective  is  equal  to  (Pew  allows  determination  of  the  run 
time  of  the  battery  at  the  chosen  power  level: 


^Effective\t=r  ~  QcUT • 


For  sizing  purposes,  the  battery’s  specific  energy  and  energy  density  may  be  determined 
from  integration  of  equation  (5): 

ES  =  ^,  ED  =  (6) 

m  v 

Through  inclusion  in  a  computer  program,  sizing  simulations  may  be  carried  out  to 
determine  available  energy  for  the  particular  power  level  demanded,  rather  than  relying  on  a 
more  approximate  value.  An  additional  benefit  of  modeling  battery  behavior  in  this  way  is  the 
handicapping  of  maximum  power.  The  use  of  constant  values  for  energy  content  requires  a 
similarly  crude  estimation  of  maximum  power  output  for  a  battery  system,  i.e.,  values  for  the 
battery’s  specific  energy,  specific  power,  and  energy  density.  For  the  model  presented  here,  the 
maximum  power  is  automatically  limited  to 

p  _  (EFULl)2  (n\ 

• Maximum  —  4  ‘  ' '  / 

The  derivation  of  equation  (7)  is  given  in  appendix  B. 

Appendix  C  presents  an  example  to  demonstrate  the  process  by  which  this  model  might  be 
used  to  predict  battery  performance. 
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APPENDIX  A 

DERIVATION  OF  EQUATION  (5) 


The  voltage  across  the  terminals  of  a  battery  is  equivalent  to  the  open-circuit  voltage  of  the 
cell,  less  the  voltage  drop  as  calculated  by  Ohm’s  law: 


E  =  E0C-(RC  •  /).  (A-l) 

Drawing  on  the  model  developed  in  Tremblay  et  al.,1  we  may  compute  the  open-circuit  potential 
as  a  function  of  the  state-of-charge  as  in  equation  (4).  Substituting  (4)  into  (A-l)  and  applying 
equation  (2)  yields 


I  ~  ^EffectiveOo6  *))P  < 

E  =  -  (SS)  +  M  •  e*P(-S  •  C»  - 


fic 


(wn)'' 


(A-2) 


Also,  as  the  original  model  presented  in  Tremblay  et  al.1  is  for  a  constant-capacity  system  with 
no  Peukert  effect,  the  tracked  capacity  in  equation  (A-2)  should  be  the  effective  capacity.  The 
effective  capacity  discharged  at  a  particular  current  (equation  (3)),  evaluated  for  a  current  that  is 
a  function  of  time,  may  be  written  alternatively  as  an  integral: 


C Effective  ~  Effective  '  At)  —  f  I  Eff  ectivei^d-t . 


If  the  current  function  is  assumed  to  be  smooth,  then  it  may  be  rewritten  as  the  derivative  of 
capacity  discharged  with  respect  to  time: 


Effective)  _  / { t }  .  (A-3) 

Removing  references  to  current  and  replacing  them  with  capacity,  we  obtain  a  final  expression 
for  potential: 


E  = 


K ’ QCUT 


Qcut  -  c Effective 


j  T  *  £xp(  B  ’  CEffective)) 


Rr 


^R£Ef£ective)  pe - 1  Pe 


(A-4) 


For  a  constant-power  system, 

Eoc  =  E„  -  +  (.A  ■  exp  (-B  ■  C))  , 


A-l 


P  =  E-l  = 


Rc  *  Roc  — 


(^Effective )  ( lpe-l'\\Pe 


or1))? 


/  ^(Cg/Tectfpe) 

dt 


P  =  E, 


oc 


Effective )  f  f pe-l^P* 


ape-l\Ve  n  (d(c Effective)  f  ,pe-l\\Pe 

0  JJ  “  ^  5  Vo  JJ 


Z,  i 

a  .  ^(CEffetlve))^  _  b  .  ^(.CEffetive)ye  +  p_Q 


where 


(z(pe-i)\ 

A)  PC  )> 


b  = 


{QCUT -effective)  +  ^  '  ®XP(  B  *  ^  ) 


(A-5) 


The  derived  function  for  a  constant-power  system  (equation  (A-5))  has  the  form  of  a 
quadratic  equation,  which  may  be  solved  for  the  value  of  the  time-derivative  of  the  effective 
capacity  raised  to  the  power  of  the  reciprocal  of  the  Peukert  effect  exponent.  As  a  quadratic 
equation,  there  are  two  solutions:  the  positive  and  the  negative  radical  in  the  numerator.  By 
inspection,  the  physically  correct  solution  to  the  differential  equation  is  the  negative  of  the 
radical  so  as  to  have  zero  current  in  the  zero-power  case. 

In  this  form,  the  equation  is  an  ordinary  differential  equation  (ODE)  and  may  be  integrated 
to  solve  for  the  effective  discharge  against  time,  assuming  an  initial  condition  for  battery 
discharge,  e.g.,  that  at  the  start  (t  =  0)  the  battery  has  not  been  discharged  (C Effective  -  0).  The 
result  is  equation  (5): 

d(Cg/y<;ctft>e)  _  fft  -  Vb2  -  4  •  a  •  P  PG 

dt  ~  [  2  a  J  ' 

R E f f ective\t-0  —  0  • 


A-2 


APPENDIX  B 

DERIV  ATION  OF  EQUATION  (7) 

The  power  delivered  by  a  cell  is  equivalent  to  the  product  of  the  potential  and  current: 

P  =  E  ■  I . 

The  cell  potential  may  be  written  as  the  open-circuit  potential  less  the  ohmic  loss  due  to  internal 
resistance: 

E  =  Eoc  —  I  •  Rc . 

The  open-circuit  voltage  is  maximized  for  a  fully-charged  cell  at  Vfuu,  making  the  equation  for 
power  a  function  of  two  constants  (open-circuit  potential  and  cell  resistance)  and  one  variable 
(current): 

P  =  1  -  (Eoc  —  /  •  Rc), 

P  =  E0C  •/  -Rc  ■  I2.  (B-l) 

Taking  the  derivative  and  setting  to  zero,  we  find  the  location  of  the  maximum: 

%  =  £oc-2  -Hc -I, 

0  =  Eoc  —  2  •  Rc  •  I , 
i  —  goc 

‘Maximum  Power  2-  Rc' 

Substituting  into  equation  (B-l),  we  determine  the  maximum  power  and  we  recover  equation  (7): 

P  —P  gOC  n  (  EOC  'i2  _  (gfULL)2 

^Maximum  -  *OC  '  2.Rc  « C’\2.Rc)  ~  4  R(.  * 


B-l  (B-2  blank) 


APPENDIX  C 

APPLICATION  OF  MODEL  TO  AN  ACTUAL  BATTERY 


To  demonstrate  the  process  by  which  this  model  might  be  used  to  predict  battery 
performance,  the  following  example  is  presented.  To  begin,  a  battery  cell  of  interest  is  selected 
and  its  data  sheet  is  analyzed  to  determine  the  parameters  defined  in  table  1  in  the  main  text.  It  is 
important  to  state  that  determination  of  the  battery’s  parameters  is  a  subjective  and  inexact 
process.  The  same  data  sheet  may  produce  two  slightly  different  sets  of  values  when  analyzed 
by  two  different  engineers  as  this  battery  model  is  an  idealization  of  the  real  discharge  behavior 
of  battery  cells. 

For  this  example,  the  chosen  cell  is  the  Saft  Li-Ion  VL  52  E  cell  (see  data  sheet  at  the  end  of 
this  appendix).  The  values  of  the  parameters  for  this  cell  are  listed  in  table  C-l. 

Table  C-l.  Safi  VL  52  E  Mode l  Parameters 


Symbol 

Battery  Parameter 

Efull 

4.1  volts 

Eexp 

3.9  volts 

Exom 

3.2  volts 

Ecut 

2.5  volts 

h 

48.9  amperes 

m 

1.0  kg 

pe 

1.035 

Qexp 

2.5  A-hr 

Qnom 

45  A-hr 

Qcut 

48.9  A-hr 

Rc 

2e-3  ohms 

V 

0.48  liter 

The  mass  and  volume  of  the  cell  are  explicitly  defined  in  the  data  sheet  as  1 .0  kg  and  0.48 
liter,  respectively.  By  inspection  of  the  C-rate  discharge  curve*  at  25°C  on  the  back  side  of  the 
data  sheet,  the  value  of  Efull  is  4. 1  volts,  Eexp  is  4.0  volts,  Esom  is  3.2  volts,  and  Ecur  is  2.5 
volts.  The  corresponding  capacities  are  approximately  Qexp  ~  2.5  A-hr,  Q\om~  45  A-hr,  and 
Qcirr=  48.9  A-hr.  As  the  discharge  is  listed  as  the  C-rate,  the  value  of  /<>  is  equivalent  to  Qcur 


*“C-rate”  is  an  industry  term  for  the  current  that  will  completely  discharge  the  cell  in  1  hour.  Other  discharge 
currents  arc  expressed  relative  to  this  value;  c.g.,  “C/10”  is  one-tenth  of  the  C-rate,  or  the  rate  at  which  the  cell 
should  fully  discharge  in  10  hours,  barring  any  effect  of  discharge  rate  on  capacity.  These  values  arc,  obviously  , 
approximations. 


C-l 


Table  C-2.  Saft  VL  52  E  Discharge  Rale  and  Capacity  Calculated  for  pe  =  1.035 


Discharge  Rate 
(C-Rate) 

Data  Sheet  Capacity 
(A-hr) 

Calculated  Capacity 
(A-hr) 

C 

48.9 

48.8 

C/2 

50.0 

50.0 

C/3 

50.9 

50.7 

C/5 

51.7 

51.6 

C/7 

52.0 

52.2 

C/10 

52.9 

52.9 

divided  by  1  hour,  48.9  A.  The  value  of  the  Peukert  effect  exponent  is  estimated  by  measuring 
the  cell  capacity  against  the  discharge  rate  from  the  discharge  curves  on  the  data  sheet.  The 
values  measured  from  the  data  sheet  and  the  calculated  values  for  the  estimated  Peukert  effect 
exponent  of  1 .035  are  shown  in  table  C-2.  Lastly,  the  value  of  the  internal  resistance  of  the  cell 
Rc  is  estimated  from  the  observed  initial  drop  in  the  C-rate  discharge  curve.  (In  many  cases  a 
value  of  internal  resistance  is  provided  by  the  manufacturer  and  estimation  is  not  required.)  At 
zero  discharged  capacity,  the  cell  voltage  decreases  from  4.1  to  4.0  volts.  By  Ohm’s  law,  the 
difference  in  voltage  is  equal  to  the  product  of  the  current  and  the  resistance,  so  Rc~  0.1  volt/ 

/o  ~  0.002  ohm. 

To  solve  equation  (5),  a  short  program  was  written  for  MATLAB  that  employed  the  solver 
ODE45.  To  avoid  integrating  over  an  uncertain  time  interval,  the  ODE  is  inverted  and  integrated 
over  the  capacity  interval  of  0  to  Qcut-  Derivation  of  this  format  is  as  follows: 


d (^Effective )  _ 

dt 


b-\b2-4-a-P 


2  a 


pe 


dt 


d(C Effective ) 


dt 


_  T _ 2  a _ _  T _ 2  •  a _  b  +  Vb2  -4  a  P|pf 

lb  -  Vb2  -  4  •  a  •  pj  [b  -  Vb2-  4  ■  a-  P  b  +  Vb2  -  4 • a  - PJ 
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The  time  calculated  as  corresponding  to  the  value  of  Qcur  allows  for  the  determination  of 
the  cell  energy  as  a  function  of  power.  While  this  is  a  useful  tool  for  estimating  battery 
performance,  a  conservative  approach  requires  imposing  additional  limits  on  the  battery 
operation. 

The  cell  used  in  this  example — the  Saft  VL  52  E — has  a  stated  maximum  discharge  current 
of  52  amperes,  corresponding  approximately  to  the  1  -C  rate.  The  published  specific  energy  is 
also  given  as  185  W-h/kg.  Figure  C-l  plots  the  results  of  this  model  for  the  VL  52  E  both  with 
and  without  the  power  and  energy  limits  on  the  cell  specified  by  the  manufacturer. 
Programatically,  these  bounds  may  be  applied  using  a  pair  of  “IF”  statements;  i.e.  if  the  required 
power  exceeds  the  limit,  the  energy  is  zero,  and  if  the  calculated  energy  exceeds  the  limit,  it  is 
reduced  to  the  limiting  value. 


Figure  C-l.  VL  52  E  Modeled  Performance  With  and  Without  Manufacturer  Limits 
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Rechargeable  lithium-ion  battery 
VL  52  E  -  high  energy  cell 


•  Excellent  power  density  ond 
specific  energy 

•  100%  ooubmbc  efficiency 

•  Low  mointenance  battery 

•  Long  cycle  fife 

•  No  memory  effect 

Typioll  ipplicitiona 

•  High  energy  applications 

•  Defense 

•  Space 

Key  feturee 

•  Graphite-based  anode 

•  Nickel  alloy  oxide-based  cathode 

•  Sold  only  as  assembled  batteries 

•  Incorporation  of  electronics  for 
performance  efficiency 

-  Cherge/floabng/discherge 
management 

-  Cen  balancing 


Battery  level  safety 
•  Incorporation  of  several  levels  of 
redundant  safety  features  to 
prevent  abuse  conditions  such  as 
overcharge,  over  discharge  and 
short  circuit 


Electrical  cfraract»rtTtic» 


Nominal  vofcege 

38V 

Lower  voltage  h me  for  discharge 

25V 

Mmmum  discharge  current  at  RT  ooronuous 

52  Amps 

Nominal  capacity  at  4  1  V/2  5  V  and  25*  C 

52  [C/7)  Ah 

Specie  energy- 

165  Wh/Kfl 

Energy  density- 

385  Wh/I 

Mechanical  charact  ana  tics 

Nominal  diameter 

54  mm 

Nominal  height 

206  mm 

Norrnnal  weight 

1  Okg 

Nominal  volume 

- U"2BT 

Call  operating  condition* 

Charge  method  Constant  current/conscant  voltage  (CCCV) 

Charging  voltage 

4  1*004V 

Recommended  continuous  charge  current  at  25*  C 

C/7 

Operating  temperature  Charge 

Discharge 

♦  5*  C  to  +35*  C 
-  25*  Cto*  55*  C 

Storage  and  transportation  temperature 

-  40*  C  to  ♦  85*  C 

End  of  charge  detection 

150  mA 

Total  RT  charging  tine 

10  hours 

hduo&c  termnatc 

'Charge  to  4. 1  V  ana  C/IQ  rats 

August  2008 
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Rats  capability  at  25'  C  from  4.1  V  to  2  5  V 
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Rato  capability  at  C/5  from  4.1  V  to  2.5  V 


Capacity  (Ah) 
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